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Abstract 

We have reanalyzed the it^p scattering data at low energy in the Coulomb- 
nuclear interference region as measured by the CHAOS group at TRIUMF 
with the aim to determine the pion-nucleon a term. The resulting value 
a = (43 ± 12)MeV is significantly lower than other recently published values. 
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1. Introduction 

In 2008 Ellis, Olive and Savage published a paper [1] on the elastic scat- 
tering of supersymmetric cold dark matter particles on nucleons and pointed 
out that the cross sections depend strongly on the value of the pion-nucleon 
sigma term cr^-jy. This is but one example of the role of the sigma term, a 
concept that was introduced in chiral perturbation theory (ChPT) to mea- 
sure the explicit breaking of chiral symmetry due to non-zero masses of light 
quarks [2]. The sigma term represents the contribution from the finite quark 
masses to the mass of the proton. Its value is related to the strange quark 
content of the nucleon. The pion-nucleon sigma term a^N is also related to 
the value of the pion-nucleon invariant amplitude at the unphysical Cheng- 
Dashen point where s — u = 0, t = 2m^. (Here, s, t, u are the Mandelstam 
variables). Consequently, its determination is mostly attempted through 
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pion-nucleon scattering experiments at low energies and by spectroscopy of 
the hadronic shift and width in the ground state of pionic hydrogen [3]. 

Ellis et al. [1] pleaded strongly for an experimental campaign to better 
determine the pion nucleon sigma term. Indeed, due to the difficulty of 
experiments with low energy pions, which are notoriously plagued by pion 
decay and the emerging muon background, the resulting cross sections from 
different tcN scattering experiments often do not agree with each other and 
the derived values of ct^n differ substantially, often by more than the quoted 
systematic errors. Ironically, however, the experimental campaign requested 
by the authors of [1] already existed and had been published [4] after a strong 
and long running effort with the CHAOS detector at TRIUMF [5], which was 
a dedicated detector system developed to cope with the peculiarities of low- 
energy pion scattering. The experiment provided differential cross sections 
for elastic scattering of positively and negatively charged pions off hydrogen 
at five energies between 19.9 and 43.3MeV in fine angular steps ranging from 
the Coulomb-nuclear interference region to nearly 180°. 

What was really missing in ref.[4] was a theoretical analysis leading to 
a value of cy n ^. The authors of ref. [4] merely considered the isospin-even 
forward scattering amplitudes ReD + (T n ) obtained directly from 7r + p and n~p 
differential cross sections and extrapolated the values to threshold (T n = 0) 
using the functional forms given by the KH80 phase shifts [6] or, alternatively, 
from the more recent SAID FA02 phase shifts [7]. In both cases the threshold 
values and the related isospin-even scattering lengths turned out to be smaller 
than in the previous analyses. From this observation it was qualitatively 
concluded that the CHAOS data favour values of a^N that are smaller than 
recently claimed [8]. 

It is the purpose of this paper to pursue this line of reasoning and to 
derive the value of the irN sigma term using the CHAOS cross sections 
as much as possible. This approach differs substantially from the method 
favored by the GWU/TRIUMF group [8] who uses the huge nN data base 
of SAID hoping that the errors of partially contradictory measurements are 
averaged out. We, instead, prefer to rely as much as possible on the results 
taken from the most advanced pion spectrometer to date. Furthermore, by 
adopting the analysis methods and the notation of the Karlsruhe group [11], 
we use constraints that warrant analyticity and unitarity. For a criticism of 
the VPI/GWU methods see [12]. 
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2. Formalism 



The ttN a term is defined as a matrix element of the quark mass term: 

m u + m d j 

a = \p\uu + dd\p), 

4m 

where m u and are up and down quark masses, and m is the proton mass. 
It is an empirical measure of the chiral asymmetry generated by the up and 
down quark masses. The a term may be written in the form: 

m u + rrid (p\uu + dd — 2ss\p) 

a = , 

Am 1 — y 

where the parameter y, the strange quark content of the proton, is defined 
as: 

(p\ss\p) 

y - 



(p\uu + dd\p) 



In a one loop calculation of chiral perturbation theory (ChPT) Gasser et al. 
[9] obtained: 

35.5±5 w Ar 
a = MeV. 

1-1/ 

This relation provides a simple way to calculate the strange quark content 
of the proton from the known value of the a term. The low energy theorem 
[10], derived within perturbation chiral symmetry, relates the so called ex- 
perimental 7iN sigma term E to the isospin-even scattering amplitude D + at 
the Cheng-Dashen point v = 0, t = 2m\: 



^ = F^D + (u = 0,t = 2m 2 n ). 

Here F% = 92.4 MeV 2 is the pion decay constant. We use the notation from 
ref. [11]. The pion mass is denoted by m^, the nucleon mass by m, and s, 
u, t are Mandelstam variables and v — (s — uj/Am^ The amplitude D + is 
defined in terms of the ttN invariant amplitudes A and B: D + = A + + vB + . 
D + denotes the D + amplitude from which the pseudo vector Born term is 
subtracted: D + = D + — D^ pu . Within the framework of ChPT Gasser et al. 
[2] obtained: 

cr = S - A a - A R , 
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A CT = (15.2 ± 0.4)MeV, A R = 0.35MeV. 

Since the Cheng-Dashen point is outside the physical region, one has to 
perform an analytic continuation of the D + amplitude to this point. Man- 
delstam analyticity, unitarity and crossing symmetry of the ttN invariant 
amplitudes are very strong constraints when analyzing experimental data or 
when performing an analytic continuation of the invariant amplitudes out- 
side the physical region. The most frequently used method for that purpose 
is the application of dispersion relations along different curves in the Man- 
delstam plane. It is important to stress that use of dispersion relations as 
a method of analytic continuation outside the physical region assumes input 
from the whole energy region in the physical irN channel. Using results from 
phase-shift analysis (PWA) in dispersion relations implies the use of a whole 
body of pion-nucleon data. Two different methods, both based on dispersion 
relations, are used for calculation of the D + amplitude at the Cheng-Dashen 
point. In a first approach, the dispersion curves passes through the Cheng- 
Dashen point which allows calculating the amplitude D directly [14, 16, 17]. 
In a second approach [8, 13], one determines the coefficients in the Taylor 
expansion of D + around the center of the Mandelstam triangle (also known 
as a sub-threshold expansion). Invariant amplitudes are real inside the Man- 
delstam triangle. The D + amplitude is crossing symmetric and is a function 
ofz/ 2 : 

m,n 

At the Cheng-Dashen point one has: 

S = F^/3 + (0,2m 2 ) = F^(4 + 4 1 •2m2 + ...) 
E = E d + A* 

Erf stands for a leading contribution, linear in t. The term Ad, a curvature 
term [2], describes contributions of higher order in t. Calculations show 
[2, 16] that the curvature term is determined mainly by contributions from 
the t-channel, and is considered as a known quantity. In our calculation 
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we use the value from [2] A^ = 11.9 MeV. The coefficients cfgo and dgi are 
given in terms of the D + amplitude and its forward slope at the center of 
the Mandelstam triangle: 

4 = D + (0,0),4 1= ^M . 

01 t=o 

It is worthwhile to point out that for t = 0: C + (v,t = 0) = D + (u,t = 0) 
and (Jq = Cq , d^i = Cq~i- The amplitude C + is mentioned because it is used 
in our approach to extract information from the low energy data. Values 
of the a term [8, 13, 14, 15, 16, 17] range from about 50 MeV to about 75 
MeV, corresponding to large values of the strange quark content y ranging 
from 0.3 to 0.5. The results depend on the method and techniques used to 
extrapolate to the Cheng-Dashen point and on the input used. 



3. Method 

Whichever of the methods mentioned above is used, data from low energy 
7i^~p scattering assure more stable and more reliable extrapolations of the C + 
amplitude to the Cheng-Dashen point. The goal of the CHAOS [4] experi- 
ment was to obtain high quality data needed for that purpose. In a CHAOS 
experiment 7r ± p differential cross sections were measured at low energies and 
at small angles in the Coulomb-nuclear interference region, where the known 
Coulomb non-spin-flip amplitude G c [6, 18] interferes with the nuclear am- 
plitude G + . The real part of forward nuclear amplitude C + (D + ) is obtained 
from the experimental scattering data using the formula: 

ReC + {q 2 ,t = 0) = ^^-ReG + (q 2 ,t = 0) = 
m 



lim 

*->o m 



da 4. da 



dQ. dQ. 

4ReG c (t) 



limA + (t), 



where Gc is the known Coulomb non-spin-flip amplitude. Measurements of 
differential cross sections for 7r^p scattering in the Coulomb-nuclear interfer- 
ence region allow a determination of the real part of the forward amplitude 
C + at five energies covered by the CHAOS experiment in the low energy re- 
gion ranging from T n = 19.9 MeV to 43.3 MeV. In a first step of our analysis 
the values A + , derived from the CHAOS data at a given energy, are fitted 
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to a polynomial in t and extrapolated to the forward point t = 0. A robust 
convergence test function method [19] is used in which an additional term $ 
is added to x 2 i n order to assure a soft cut off of higher terms in a polynomial 
expansion: 

$ is of the form $ = A Y2^=o a n( n + l) 3 ; where A is a smoothing parameter. 
Results for ReD + are shown in Fig.l. 
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Figure 1: Data points: Real parts of the forward D + amplitude from the CHAOS 
experiment [4]. Solid curve: Fit curve as described in the text. 

If the real part of the irN amplitude C± is consistent with the forward 
dispersion relations (i.e. fixed t analyticity) then, according to so called 
Schnitzer-Salzman plot [20, 21], the quantity: 

Y±(u) = ReC ± (u) - I±(w) - C N ±(u) 

is linear in u, energy of the incoming pion in the lab. system. In the above 
formula Cn ± is a Born term and I± stands for an integral: 

oo 



doj'. 
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Figure 2: Schnitzer-Salzman plot for the CHAOS data (see text). 



Figure 2 shows a quite satisfactory behavior of the function Y obtained from 
the CHAOS data. In the next step, the obtained values of ReC + (q 2 ,0) are 
fitted to a polynomial of second order in q 2 : 

ReC + (q 2 ,0) = c + Cl q 2 + c 2 q\ (1) 

and extrapolated to the threshold (q 2 — 0, t — 0). The first two coefficients c$ 
and Ci determine ReC + (q 2 ,t = 0) and the derivative aC+ ^ q 2~^ at threshold, 
respectively. Using the partial wave decomposition of ReC + and the effective 
range formula for partial waves fi± = Re [j^j — <? 2i (<^± + h±q 2 ) close to 

the threshold, one obtains the following approximation of the C + amplitude 
[11]: 

Am 2 ) 4vr(l + x) 0+ 2\ l+ x ~ Am 2 ) 

+ (ft + 2< ++ < + ^_).^ 

+ higher order terms in t and q 2 . 

Here, isoscalar s-and p-wave scattering lengths, is the s-wave 

effective range parameter [11] and x = —. 
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The second term, according to Geffens sum rule [28], is proportional to a 
first derivative of the C + amplitude with respect to q 2 at threshold: 



1 dC + (q 2 ,0) 



4tt(1 + x) dq 2 



&+_ + 2a+ + + a+_ + 



q 2 =0 



2m ■ ra„ 



2af + + a\_ 



Cl 



47r(l + x) 



2m • 7T7.Tr 



(2) 



Taking a derivative with respect to t (forward slope), one obtains: 
1 8C + {Q,t) 



47r(l + x) <% 



= \ + + + «+_)■ (3) 



From formulas (2) and (3) one may derive an expression relating 



dC+(q 2 ,t=0) 
dq 2 



q 2 =0 

dReC + (v = m n , t) 



to the forward slope 



dC+(g 2 =0,t) 
dt 



t=0 



dt 



Cl 



t=0 



2 V ; V 0+ 2m-m, 



The corresponding forward slope of the amplitude C + is obtained by sub- 
tracting a pseudovector Born term: 



dReC + (v = mm t) 



dt t=o 
dReC + (v = m w , t) 



dt 



dReC% pv (v = m w ,t) 



dt 



t=0 



The Born term and its derivative are explicitly known [11]: 



dReCjj pv {v = m^t) 



dt 



9 w 



4m 3 uj" — uj 



t=o 



2 _ , ,2 



uj 



UJ + UJB 



In order to determine the coefficient c?qi = c 01 



+ _-+_ ( dC+(v=0,t) 

at 



and subtracts the forward slope at threshold: 



one adds 



t=o 



dReC + (v = 0,t) 



t=o 



dReC + (v = Q,t) 
9t 



dReC + (v = mis, t) 



dt 

dReC + (u = mis, t) 



+ 



dReC + (u = m w , t) 



dt 



t=o 



t=0 



dt 



Diff 



t=o 



01 



+Diff 



dReC±(v = m^t) 



dt 



t=o 



01 



(4) 



Ai and Diff stand for a difference between the forward slope of C* + and C + 
respectivelly at the center of the Mandelstam triangle (u — 0, t — 0) and at 
threshold (z/ = m n , t = 0). From dispersion relations for the forward slope of 
C + [11] one obtains the following expression: 



Ai 



dReC% pu (v = m^t) 



9* 



da;' 2a;' + a; 



7T 



t=0 



da/ 



a/(a;' — m 2 ) 



2m / a/ 2 (a/ — m^) 2 



dReC£ pu (v = m n ,t) 



imC + (a/,0) 



Ot 



(5) 



Both integrals are fast converging and may be accurately evaluated using re- 
sults from existing partial wave analyses or, by virtue of the optical theorem, 
from existing data for total cross sections. 



4. Results 

The CHAOS data [4] consist of differential cross sections grouped in two 
angular regions. The first group contains data at forward lab angles, typically 
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below 35°, which were covered by a range telescope for particle identification. 
The second group comprises all larger angles. We have performed two sep- 
arate analyses of the CHAOS data. In the first analysis we have included 
data at forward angles only. Values of A + (t) obtained from the experimen- 
tal data have been fitted to a polynomial of third order in t. In the second 
analysis, all data were included and fitted to a polynomial of order five. In 
such a way two sets of values of ReC + were obtained from the CHAOS data 
and averaged. This way we deliberately increased the weight of the small 
angle data, keeping in mind the design features of the CHAOS detector. The 
resulting values of ReC + are nuclear. Hadronic values were calculated ap- 
plying electromagnetic corrections according to the Nordita procedure from 
ref. [18] and used as an input for determination of the coefficients d,Q Q and 
u 01 . 

Forward dispersion relations (FDR) are an essential part of our extra- 
polation of the C + amplitude to the threshold. In the first step, values of 
ReC + , obtained from the CHAOS data, combined with the once subtracted 
dispersion relations for the C + amplitude, are used to calculate its value at 
the threshold. The FDR for the amplitude C + , subtracted at the threshold, 
read [11]: 



where k is the pion lab momentum, f 2 is the ttN coupling constant and 
the nucleon Born term. For each energy where ReC + is determined from the 
data, one calculates the dispersion integral and obtains the corresponding 
value of the subtraction constant C + (m n ). Partial wave analyses do not 
give errors of partial waves so that the error of the integral in (6) may not 
be given and errors quoted in our results are due only to the numerical 
procedure. Our calculations show that errors of the subtraction constant 
due to the uncertainty of dispersion integrals are negligible compared to the 
experimental errors. 

As a result from the first step we obtain an average value of the sub- 



ReC + (u 



C + ( m7T ) + C+(Lo)-C+( m7T ) 

oo 

2k 2 f ImC+{u') 




(6) 



ml lo 2 — Wr ' ' 2m ' 



10 



traction constant with the corresponding error. We stress that the obtained 
average value is not our final result. It is used as a part of the input when 
calculating final values of the amplitude and its forward slope at threshold. 
Forward dispersion relations for the amplitude C + evaluated at the center 
of the Mandelstam triangle, give useful relations connecting values at the 
threshold and at the center of the Mandelstam triangle: 



OG 



7r J u'(u' 2 — ml) 

m w 

and 

^ + (0) = 4o = C + (m„) + A 2 = C + {m^) + 1.88/ 2 + A 2 . (8) 

The dispersion integrals in (5) and (7) may be accurately evaluated using 
available scattering data. In the energy region where results from PWA exist 
(k < k max ) imaginary parts are calculated from partial waves (GWU/VPI: 
k m ax = 2.6 GeV/c; Ka84: k max =6.0 GeV/c). Tables of total ir^p cross 
sections, needed to calculate ImC + , are available up to lab momenta of 
k = 340 GeV/c [22]. Parametrization of the total cross sections at high 
energies are also available [25]. The integral in (7) is fast converging. More 
than 97% of its value is due to contributions below k = 2.6 GeV/c, so that 
uncertainties from the high energy part of the integral may be neglected. 
Using the GWU/VPI solution Fa08 [SAID] one obtains a value of A 2 = 
— 1.381m" 1 . Ai has been calculated using the same input. Our evaluations, 
using several PWA solutions (Ka84, Fa02, GW06, Fa08 ), show that the 
integral l\ is saturated already at co max which corresponds to the highest lab 
momentum in the GWU/VPI partial wave solution. Using the discrepancy 
method we found that contributions from higher energies to the integral I\ do 
not exceed one percent of its value. Hence, uncertainties from high energies 
may be neglected. Using results from the GWU/VPI PWA solution Fa08 
and tables for total cross sections at high energies, a value of Ai = -133.8 
GeV -3 = -0.364 ml was obtained. Extrapolation to threshold yields the first 
two coefficients in expansion (1): 

c = (-0.140 ±0.013)m; 1 
c t = (2.146 ±0.187)m; 3 , 

which are the main result of our analysis. 
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The only parameter left to be determined is the s- wave effective range 
parameter 6q + . One may follow the authors of ref. [23] and use the Karlsruhe 
value &o+ = (—0.044 ± 0.007)m~ 3 from ref. [11]. The very same value was 
obtained from partial wave relations derived from the fixed t dispersion rela- 
tions [24]. Due to the fact that partial waves from partial wave relations are 
strictly consistent with analyticity, the method allows a reliable determina- 
tion of threshold parameters. Unfortunately, like some other programs from 
the Karlsruhe group, Koch's program for evaluation of partial wave relations 
was lost. To our knowledge there are no recent determinations of the s-wave 
effective range parameters of such a degree of sophistication. 

Until results from a new evaluation of partial wave relations become avail- 
able, we rely on a simple extrapolation of the s-wave to threshold. The value 
we use in the present analysis is an average of the Karlsruhe value and a 
value obtained from the current GWU/VPI partial wave solution: 



Inserting values of A 1; A 2 and into formulas (2), (4), and (8) we obtain: 
c+ = (-1.3780 ± 0M30)m-\ 2a+ + + a+_ = (0.1994 ± 0.0136)m^ 3 , c£ = 
(1.0710 ± 0.0980)m~ 3 , and S d = (46.7 ± 12.0) MeV. Adding the curvature 
term, we obtain our result for the ttN cr-term as derived from the CHAOS 



As stated before, the quoted errors are mainly due to errors having their 
origin in the experimental data. Errors due to uncertainties of dispersion 
integrals, as shown by calculations, are negligible compared to the errors 
having their origin in the experimental data. 

5. Discussion 

The result obtained for the coefficient Cq corresponds to a s-wave scat- 
tering length ciq + = (—0.0097 ± 0.0009)m^ 1 which is comparable to the old 
Karlsruhe result [11]. The value of combination of p- waves 2ai + +a{_ has not 
changed significantly during the last three decades. Our value 2af + + a+_ = 
(0.1994 ± 0.0136)m~ 3 is in good agreement with those obtained in partial 
wave analyses performed during the last several years [26, 27]. It is coupled 
with the s-wave effective range parameter b$ + through Geffen's sum rule [28]. 




(-0.050 ±0.004)m; 3 . 



data: 



S = (58.6 ± 12.0) MeV, a 



(43.1 ± 12.0) MeV. 
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To study the dependence of our results on the s-wave scattering length a^ + , 
we alternatively used the newest result from the pionic hydrogen [3] exper- 
iment as part of the input in our fit, taking a$ + = (0.0044 ± 0.0022)m~ 1 . 
Although the resulting fit to the data is numerically quite satisfactory, the 
obtained values are in strong disagreement with other low energy parame- 
ters. Namely, the value of C\ = (1.31 ± 0.18)m~ 3 obtained in such a way is 
in disagreement with recently obtained values of 2af + + af_ unless a large 
negative value b$ + — 0.09m~ 3 for the s-wave effective range is allowed. 

6. Conclusions 

Our value of the irN <r-term, based on an analysis of the CHAOS data, 
is significantly lower than most values published in the last ten years. Given 
the systematic uncertainties of the CHAOS experiment one could not expect 
a precision determination of its value. Using an analysis which respects the 
analytic properties of the ttN amplitudes we were, however, able to confirm 
quantitatively the conjecture of Denz et al. [4] that the sigma term is small, 
comparable to the canonical value a = 49 MeV of Koch et al. [6, 14]. This 
is definitely much smaller than the result of the phenomenological extraction 
by the SAID group [8]. With no experimental facility available to improve 
the ttN data base in the forseeable future it was important to question that 
result. Among other consequences we point out that, using our result, the 
difficulty with the unreasonably large strange quark content of proton no 
longer exists. 
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